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Abstract. The fluctuations of macroscopic observables in quantum systems 
which are in a nonequilibrium steady state are studied rigorously in the ther- 
modynamic limit. In particular, the nonequilibrium steady state (NESS) of 
a quantum spin system that is composed of two infinite reservoirs at dif- 
ferent temperatures coupled in a bounded region is considered. Under suit- 
able assumptions on the interaction and the asymptotics of the fluctuations, a 
nonequilibrium central limit theorem for the fluctuations of the NESS expec- 
tation value of the empirical average of observables is proven by showing the 
analyticity of the corresponding moment generating function in a neighbour- 
hood of the origin using the cluster expansion method. Furthermore, the large 
deviation principle is proven for a class of observables by establishing the ex- 
istence and differentiability of the corresponding moment generating function. 



Dedicated to Jiirg Frohlich on the occasion of his sixtieth birthday. 



1. Introduction 

The central limit theorem and large deviation theory have a long history in 
equilibrium statistical mechanics, where they play an essential role in under- 
standing the fluctuations of macroscopic observables in systems at equilibrium 
(;see for example [T] ). For instance, it is well-known that the relative entropy den- 
sity governs the large deviations of empirical measures for Gibbs random fields in 
classical statistical mechanics, see for example [T] and [2]. Noncommutative cen- 
tral limit theorems have also been established for the fluctuations of macroscopic 
observables in equilibrium quantum statistical mechanics. [3j HJ El E] Recently, the 
large deviation principle has been proven for fluctuations of some observables of 
quantum systems at equilibrium. For continuous quantum systems, the large de- 
viation principle (and a central limit theorem) for the fluctuations of the particle 
density has been established for ideal quantum gases in [7], and dilute quantum 
gases in [8], using the cluster expansion method. Furthermore, the large deviation 
principle and a central limit theorem has been shown to hold for fluctuations of 
onsite observables in quantum spin systems at equilibrium in [5] using a high tem- 
perature cluster expansion, while (part of) the large deviation principle has been 
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proven for a general class of observables of quantum spin systems at equilibrium 
in [10] using the Gibbs condition. 

Unlike in systems at equilibrium, fluctuations in nonequilibrium steady states 
(NESS) is poorly understood. Recently, there has been encouraging progress in 
understanding fluctuations of observables in specific classical statistical systems 
at NESS in the hydrodynamic limit, particularly in simple exclusion processes, 
[TT| [12] . and in stochastic Hamiltonian systems, [TH [TJ], based on ideas from 
large deviation theory. There is no analogue of the latter results in nonequi- 
librium quantum statistical mechanics starting from a (time-reversal invariant) 
microscopic description of quantum systems. □ 

In this note, we consider the nonequilibrium steady state (NESS) of a simple, 
yet nontrivial, quantum system, a quantum spin system consisting of two infi- 
nite reservoirs at different temperatures which are coupled in a bounded region, 
and we study in a mathematically rigorous manner the fluctuations of the NESS 
expectation value of the empirical average of observables in the thermodynamic 
limit. In particular, under suitable assumptions on the decay of interactions and 
the asymptotic evolution of the fluctuations, we prove the existence and ana- 
lyticity of the corresponding moment generating function in a neighbourhood of 
the origin. By Proposition 1 in [16], this implies that the fluctuations satisfy a 
nonequilibrium central limit theorem. In the special case when there is no cou- 
pling between the reservoirs, the state of each reservoir is an equilibrium KMS 
state, and we recover a central limit theorem for the fluctuations of general ob- 
servables in quantum spin systems at equilibrium, which was mentioned in [9] and 
announced in |10| . but which has not been explicitly proven so far. Furthermore, 
the differentiability of the moment generating function for a class of observables 
implies the large deviation principle by the Ellis- Gartner theorem, [2J. 

The organization of this note is as follows. In Section 2, we discuss the model 
system we consider, recall the definition of NESS, and describe the main problem. 
In Section 3, we state the main result, which we prove in section 4. Our analysis 
relies on the algebraic formulation of quantum statistical mechanics, [T71 ITS] , the 
notion of nonequilibrium steady states, [Hi [20], large deviation theory, [2], the 
Kotecky-Preiss criteria for the convergence of the cluster expansion, [25, 26J, and 
the analysis in [9l [10] of fluctuations of observables in quantum spin systems at 
equilibrium. In order to get spatial (and temporal) profiles and to make contact 
with the quantum version of the results in [UJ, [121 EES EHh one needs to take 
the hydrodynamic limit of the reservoirs rather than the thermodynamic limit. 
Further physical applications and a discussion of the hydrodynamic limit are the 



fluctuations for stochastic quantum systems out of equilibrium has been studied in [15] 
using a phcnomcnological approach. 
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subject of future work. We note that combining the results of |2T] and p], one can 
easily extend the analysis of this note to continuous fermionic quasi-free quantum 
systems in NESS. 



1.1. Notation. Consider a d-dimensional lattice Z d . For a = (oi, • • • , a^) G Z d , 
such that aj > 0,i = 1, ■ • • , d, we define A a C Z d as the parallelepiped 

A a := {x = (xi, • ■ ■ , x d ) G Z d : < x { < a h i = 1, • • ■ , d}. 

The set of translates of A a by na := (niai, • • • , n d ad), G Z, ? = 1, • • • , d, form 
a partition P a of Z d . For A C Z d , we let N^(a) be the number of sets of this 
partition which have a nonempty intersection with A, and we let N^(a) be the 
number of sets of this partition which are included in A. Finite subsets A C Z d 
are defined to tend to Z d in the sense of van Hove, A /* Z d , if 

lim ^w_, 



for all partitions P a -[IZl EH] 

2. The Model and Description of the Problem 

2.1. The Model. We consider a quantum spin system composed of two infinite 
reservoirs, Rj, i — 1, 2, at different temperatures which are coupled in a bounded 
region AoO For the sake of concreteness, we assume that each reservoir, Ri and 
R2, occupies Z d . 

We associate with every x € Z d two finite dimensional complex Hilbert spaces 
Hi and H^. Furthermore, for X C Z ' we define the local Hilbert spaces 

We denote by U l X) i = 1,2, the kinematical algebra of local observables, the C*- 
algebra defined by U l x = B(TC X ), the algebra of bounded operators on Ti, l x ,i = 

1.2. If X fl X' — <f>, then H l XnX , = H' l x <g> H x ,,i = 1,2, and U l x is isomorphic 
to U x ® V x ,, where l l x , is the identity operator on 7i x ,,i = 1,2. Note that, by 
construction, U l x C U l x , if X C X', and [U X ,U XI ] = whenever XflX' = 0, where 
[■, ■] stands for the commutator. Furthermore, [U X ,U X ,] = for every X, X' G Z d . 
We denote by U x := IA\ ®U X . The quasi-local C*-algebras 14,14^, i = 1,2, are 
given by 

W = U X:X , cZ M x ®U x ,,U Ri = U XcZ M x ,i = 1,2, 



2 We refer the reader to [T7] for a detailed exposition of quantum spin systems. 



4 WALID K. ABOU SALEM 

where (•) denotes the norm closure. The quasi-local C*-algebras U and have 
a common identity element 1. 

In order to specify the dynamics of the uncoupled system, we introduce the 
interaction <3>* in each reservoir Rj,i = 1,2, which is defined on finite subsets of 
Z d . For A C Z d , such that |A| < oo, 

¥ : A -> W A , 

such that <3> l (A) is a selfadjoint element of W^. We define the local Hamiltonians 

H Rl (A):= y E* i (X),i = l,2. 

XcA 

We also introduce the interaction * corresponding to the local coupling of the 
reservoirs, 

* : Z d D A -> W A , 

such that *(A) is a selfadjoint element oiU\ and *(A) = if A n A = 0. The 
interaction Hamiltonian corresponding to the local coupling is formally given by 



ACZ d 

AnA ^</> 



We make the following assumptions. 
(Al) Decay of the interaction. 3A > such that 



\&\\ x ■= sup J2 e1lX E II^WII < °°, i = 1,2. 



xe& n>0 X3x 

\X\=n+l 

(A2) Decay of the coupling. 

||*|| := sup E ||*P0|| < oo. 

xcjA 
X3x 

It follows from assumption (A2) that the interaction Hamiltonian V is bound, 

< |A |||*|| < oo. 



We introduce the linear operators S and So with domain V = U Ac zdU\, such that, 
for A e U\i , A' C Z d , 

2 

»=1 xcz. d 
xnA'yt<t> 
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and 

2 



i=1 xcz d 
xnAV</> 

We also introduce, for A C Z d , the linear operators 5a and 5q,a which are given 

by 

2 

5 A (v4) :=^^[i/ Ki (A),A] + ^[y,A], 
i=i 

and 

2 

5 0)A (A) :=^[^(A),A], 
i=l 

for A € V. The following lemma follows from assumptions (Al) and (A2) and 
Theorem 6.2.4 in [T7]. 

Lemma 1. Suppose assumptions (Al) and (A2) hold. Then T> is a norm-dense 
*-subalgebra of analytic elements of the (norm) closures 5 and <5 of 6 and So 
respectively. Furthermore, 5 and So generate strongly continuous one-parameter 
groups of *- automorphisms a and cto ofli, such that 



mil net — r 

and 



lim \\a\A) -e tdA {A)\\ = 0, 

A /l£ 



lim \\aUA) -e tSo ' A (A)\\ = 0, 
for all A eW, uniformly in t GK /or £ m compacts. 

We assume that the initial state of the system, ojq, is given by 

Tr HA (e-^^ H ^A) 
u (A) := hm — — 5 — -, 

v ; A/zd Tr HA e-^i ft ^. (A) 

for A e W, where > 0. The existence of the initial state follows from assumption 
(Al); see for example [T7J [18]. We note that the initial state u> is a -i nvar i an t 
on U, and the initial state of Rj, i — 1,2, is the equilibrium KMS state at inverse 
temperature Furthermore, to>o is modular H 
We make the following additional assumptions. 



3 We recall some basic definitions from operator algebras. Consider a state a; on a C*-algebra 
U. We denote by (H u , n^, the GNS-representation of U corresponding to the state u, and by 
Muj the envelopping von Neumann algebra ir^{U)" . A state uj is called modular if its extension 
to M. u is faithful, i.e., if AQ, U = for A G M.^, then A — 0. A state ry on U is said to be 
w-normal if 3p a density matrix on 7i u such that rj(-) = Ti , t-c lo (pt: U j(-)). A state 77 on W is said 
to be w-singular if 77 > A</>, for some nonnegative A and aj-normal state 0, implies A = 0. 
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(A3) Homogeneity of the interaction. Denote by r the automorphism acting on 
U AcZ dU\ corresponding to translations. If A G U\, then t x (A) G U\ +x . We 
assume that $ is translationally invariant, i. e., for AcZ d and x G Z d , 

$(A + X ) = $(A). 
(A4) Scattering endomorphisms. The limits 

7±(A) := lim o^ 1 o a* (A) 

t— >±oo 

exist in norm for all i G W. Furthermore, 7± are norm-preserving *- 
morphisms of W. 

We now recall the notion of a nonequilibrium steady state (NESS) which was 
first introduced in [19]. We define the NESS as 



(1) uj + := lim w oa'=w o7 



t— >oo 



Note that by construction u + is an a— invariant state on U. Furthermore, since 7+ 
is a norm-preserving *-morphism and ujq is a positive functional, u + is a positive 
functional, i. e., w+(A*A) > for all AeU.^ 

Remarks. (1) Since the coupling V is bounded, it follows from a Dyson series 
expansion that 

POO rti rtm-1 

u + (A) = uj (A)+y^(z) m dh dt 2 --- dt m u {[afr{V),[---[a%{V),A]]]), 

Jo Jo Jo 



m>l 



for all A G U. This expansion is useful in computing the rates of heat transfer 
and entropy production; see for example [201 |2T| [23] . 

(2) We note that (AJ^) follows for example from the assumption that there exists 
a norm-dense subalgebra Uq C U such that the C* -dynamical systems (U,a) and 
(U,Q!o) are LMJAq)- asymptotically abelian, i.e., 



/oo 
dt\\[A,a?(B)]\\ < 00, 
-00 

/oo 
dt\\ [A, 4(B)) || < 00 
-00 



for all A, B G Uq. [19] The assumption of asymptotic abelianess implies a form 
of ergodicity of the system, and it has been verified in only few physical models, 
such as the ideal Fermi gas.\ 



4 One can show for specific models that w+ is wo-singular, and that the entropy production 
of the coupled system in NESS is strictly positive, see for example [2T1 [25] , 
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(3) One can relax assumption (A4) by weakening the definition of the NESS; 
(see for example [20j[23],). Suppose only assumptions (Al) and (A2) hold. Then 
Lemma[J\ follows. We define the NESS as the limit 



in the weak-* topology, i. e., 3 a sequence {T n } such that lim n _+oo T n = oo and 



The dual U* of U is compact in the weak-* topology, and hence the set of all 
NESS is nonempty. If assumption (A4) holds, then uj + = ujq o j + is unique. 

(4) Assumption (A3) is a technical one, and it simplifies the application of the 
cluster expansion method. 

2.2. Description of the Problem. We now describe the problem we are in- 
terested in. Given a self-adjoint element A G U , A* = A, we introduce the 
probability distribution P^, which, for O a Borel measurable subset of R, is given 



where u + is the NESS defined in ([T]) and To is the indicator function of O C R. 
Note that Pa(0) denotes the probability that the expectation value of A in state 
oj + takes values in O. It follows from ([2]) and the spectral theorem that, for a 
Borel measurable function /, 



Ja(A) 

where o~(A) C R is the spectrum of A. 

For every X G Z d , we associate A(X) G U l x , such that A(X) is selfadjoint. For 
A C Z d , we define the local observable A\ := ^^^(X). We also define 



It follows that P^ A is the probability measure associated to the empirical measure 
of A. We denote the fluctuation of A^ by 

5A A := A a -lv+(A a ). 

We ask whether 





by 
(2) 



P A (0) = u J+ (To(A)) 





(3) 
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obeys the central limit theorem, i. e., 3a 2 > such that, for all z£R, 

(4) lim u + (e izWA ) = e~ z2a2/2 , 

\yz d 

where 

a 2 = lim \A\lu + ((5A a ) 2 ). 

Another question that we address is whether the sequence 

{P^ A : A C Z d } 

satisfies the large deviation principle, i.e., whether there exists a lower semi- 
continuous convex function / : R — > R, the so called rate or entropy function, 
such that 

1 

) - 

A/Z d 



(5) lim sup — log P^ A (C) < - inf I(x) , C C R dosed 



(6) lim inf — T log P-r (O) > - inf I(x) , O C R open. 

A/Z d |A| x£0 

We now state the following assumption on the observable. 
(A5) Asymptotic evolution of the fluctuation. 

s - lim a*(A A ) = A' A + R A , 

t— >CX) 

where A' A = 4r ^2 A , cA A'(X), A' : A — > W is selfadjoint such that there 
exists A' > with 

WA'Wx = sup^e nV W A '( X )W < 00 » 

ieZd n>0 X3! 

|X|=n+l 

and there exists a finite constant C independent of A such that 

lim IAIH^a-CII =0. 

A/Z d 

Assumption (A5) is trivially satisfied if A commutes with the generator of the 
time evolution of the coupled system, i. e., the observable is a conserved quantity 
of the coupled system, in which case R A = and A' = A. We remark how this 
assumption can be relaxed after Theorem [H Section [3j 



3. Main Result 



In proving the central limit theorem and the large deviation principle, it is 
often useful to introduce the moment generating function, which physicists call 
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"free energy". We denote by F the moment generating function, which is given 
by 

(7) F(a) := lim 4r logu; + (e a|ApA ) = lim 4r hgu+(e a ^^ A{x) ). 

If F is analytic is a neighborhood of the origin, Proposition 1 in [16] imply 
the nonequilibrium central limit theorem 
C 1 (IR), then, by the Ellis-Gartner theorei 
and the rate (or entropy) function is 



. Furthermore, if F exists and is in 
the large deviation principle holds 



(8) I(x) = sup(ax — F(a)), 

aeR 

the Legendre transform of F. We note that if F exists without being differentiate, 
then the least upper bound flS} holds, but the lower bound (jSJ) fails (see (2] for 
further details). We now state the main result. 

Theorem 1. (i)Suppose assumptions (A1)-(A5) hold. Then there exists positive 
constants (3 and jl, which depend on A and X appearing in (Al) and (A5) respec- 
tively, such that, for |/3i j2 | < $, the moment generating function F defined in ([?]) 
is analytic in a for a G Mp, := {z G C : \z\ < jl}, and W\ defined in satisfies 
the nonequilibrium central limit theorem Q). 

(ii) If, in addition, A' appearing in (A5) is onsite, then, for |/3 1)2 | < $, the 
moment generating function F exists and is analytic in the strip h og 2 := {z G 

C : | ^sz | < an d the large deviation principle (TJJj and fifty hold, where the 

rate function I is given by (TJj). 

We prove both claims in the following section. 

Remarks. (1) In the particular case when the interaction between the reser- 
voirs is turned off, each reservoir is at equilibrium, and claim (i) implies 
that the fluctuations of general observables in each reservoir at equilibrium 
satisfy a central limit theorem. This is a result which has been announced 
in [H] and [TDj 7 but for which no explicit proof has been provided so far. 
On the other hand, claim (ii) is a straightforward extension of the results 
of [9] to the case at hand; see also [10] . 
(2) One may relax assumption (A5). Under the weaker assumption 

s - lim a\A A \ Aint ) = A A \ A , + R A , 

where A int and K' int are fixed bounded regions ofTL containing the interac- 
tion region A , and A A and R A satisfy the conditions of assumption (A5), 



5 see for example [5] 



10 



WALID K. ABOU SALEM 



one can show that 

F'(a)= lim — -\ — r logaj + (e a|MAl "* pA \ A -0 

A/Z<* \k\k int \ 

is analytic in a in some neighbourhood of the origin. (This follows from 
claim (i) in Theorem^ and the fact that A /* 7/ in the van Hove limit; 
see proof of (i) in the following section.) Therefore, fluctuations of the 
observable in A^ nt = Z d \A int obey a central limit theorem. Furthermore, 
claim (ii) holds in the particular case when A' is onsite, which implies that 
the large deviation principle holds for fluctuations in A c int . The latter is 
satisfied, for example, when A is onsite and commutes with the generator 
of the free evolution, 5q(A x ) = 0,x G Z d , and the interaction \1/ has a finite 
range r. In this case, A x commutes with the coupling V for x G A int = 
{x G Z d : dist(x,A ) > r}, and hence 5(A X ) = 5 (A X ) + i[V,A x ] = 0, 
A' int = A inU A = A', and R A = 0. 
(3) One can prove similar results if one has specific additional information 
about the correlations of the fluctuations, such as in systems which can 
be mapped onto quasi-free nonequilibrium system. Suppose for instance 
that there exists a constant k > independent of A C Z d , such that the 
correlation function of the fluctuation satisfy 

\u>+((8A A ) n )\ < n\ ' 



| A |n-l' 



for all n > 1. In this case, one can easily show that F is analytic in 
a for |a| < k, and hence the fluctuations W\ defined in |3|) satisfy the 
nonequilibrium central limit theorem Ol). 



4. Proof of the Main Result 

4.1. Cluster Expansion. We start by proving a useful lemma. It follows from 
assumptions (Al) and (A5) that there exist positive constants (3, a and /2, which 
depend on A and A', such that 

(9) sup ^ e 2a]A| (e AIK(A)||+^||^(A)+^(A)|| _ i) < 5 . 
For A C Z d , we define 

(10) f (A) ;= e «|A| ( ^||A'(A)|| + /3||<I.i(A)+* 2 (A)|| _ ^ 
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and £(0) = 1. We denote by S the set of finite collections of subsets of Z d . For 
S = {Ai, • • • , Afe} E S, we let \S\ — k, the cardinality of S. We define the function 

\s\ 

(11) e(S|A):=£ J2 II 

i>>1 A6{A 1 ,-,A !/ } c c5 , /i=l 

for S E S and A E S, where {-} c denotes the collection of connected subsets of 
7L d . For S E S, we define 5 := U^Aj. We have the following lemma. 

Lemma 2. Suppose assumptions (Al) and (A5) hold. Then 

(12) sup V6(S|A) < a, 

ses " 

^tz c A3ai 

AGS' 

where a appears in (TJ|) and is defined in < f77)) . 

Proof of Lemma\^ The proof is based on induction in the cardinality of S E S. 
The case |5| = 1 follows trivially from (T9T) - (TTXT) . Suppose that 

(13) sup V0(S|A)<a. 

S&.W-N A3, 

Now, for \S\ =N+1, 

N+l v 

Ee(siA) = EE E n« A ») 

As, As, v=l Ae(Ai,-,A„} c cS(i=l 

AeS AeS 

TV i/ 

= E« A >E E n« A J 

A3« s>=o Ae{Ai,— ,A„} c cS\{A}/i=i 

sE«( A >E^EIK E n« A <» 

A3i u>0 ri,— ,r„ S=l {A lr -- ,Ar s }oCS\{A} 1=1 

sE« A >E^( E «< A ')E E n« A <>r 

A3, 1^>0 A' r>0 {Aj,--- ,A r }cS\{A,A'} *=1 

AeS hnh'J=<t> {AA'Ai,- Ar}c 

|A|su P:c , eSUA} E A , 3<e , 6(5\{A}|A') 

<2^^(A)e A'gs\ { A } 

A3a; 

AeS 

Together with (TT3"j) , this implies 

(14) sup ^0(S|A)<^e(A)e |A|5 

ses.isNAr+i ABx A3x 
x ^ AGS AeS 
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However, we know from (Q and ([TO]) that 

^£(A)e |A|<i < ^ e 2 «l A l( e / i l|A'(A)||+^||* 1 (A)+* 2 (A)|| _ i) < 3. 



A3a- 



The last estimate together with (1T4|) imply that 

sup ^6(S|A)<d. 



SeS,\S\=N+l 



A3x 



+ Aes 

Claim (112p follows by induction. □ 

We now introduce further notation and definitions which are helpful in applying 
the cluster expansion method. For A C Z d , we consider B a sequence of subsets 
of A, whose union is connected, 

B = (Ax,-- - ,A fc ) c , 

where the subscript c stands for Uf =1 Aj C A connected. We denote by \B\ the 
number of elements in the sequence B,i. e., \B\ = k, and let 

B := uSA r . 

We denote by C(A) the set of all such sequences of subsets in A, 

C(A) := {B = (Ai, ■ ■ ■ , A|B|) C : A D B = uJf^A,. connected}. 
and we let 

C := U A/zd C(A). 

We introduce the measure u A on C(A), A C Z d , which, for B G C(A), is given by 

j \B\ n B 

MB) := -p m s Tr HA J2H(aA'(A r )) fl (-/^(A*) - &<I> 2 (A S )), 
aim{n A ) n=0r=1 s=n+1 

where $*, % = 1, 2, is the interaction that appears in assumption (Al), A' appears 
in assumption (A5), and = A' ((f)) = 1. Furthermore, we denote fi a := 

\im A y* Z d n\. Note that 



\^{B)\<C£\[\\aA'(A 



\\^ l (A s ) + ^ 2 (A s 

n=0 r=l s=n+l 



We also introduce a real measurable function x on C(A) x C(A) which, for Bi and 
Sj G C(A), is given by 



(15) X(B 1 ,B J ) 



0, B 1 f\B l = (f) 
-1, BiHB^ 
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The motivation for these choices will become apparent in Subsection 14.21 For 
A C Z d , we introduce the partition function 

(16) ^a:=1 + E^[ E UtiW II (1 + X(^, £,•)). 

n>l ' Bi,-,B„eC(A) i=l l<i<j<n 

We let G n be the set of all graphs with n vertices, and C n C G n the set of 
connected graphs with n vertices. Given a graph G, we denote by E(G) the set 
of edges of G and by V(G) the set of vertices of G. We introduce the combinatorial 
function ip of finite sequences (Bi, ■ ■ ■ , B n ) of C(A), A C Z d , which is given by 



(17) ifiB^--- ,B nj 



1,71= 1 

i E G ec„ n(ij) e E(G) x( S *> n ^ 2 
where the product is over all edges of G. If the graph with n vertices and an 
edge between i and j whenever x(B{,Bj) ^ is connected, then the sequence 
(-Bi, •-• , £> n ) is connected. The cluster expansion allows for an expansion of 
m logZ^ in the limit A y 7L d \ see for example [24]. We will use the Kotecky- 
Preiss criteria to prove the convergence of the cluster expansion; see [251 [26] . 

Proposition 1. Suppose assumptions (A1)-(A5) hold. Then, for \a\ < ft and 

\Pi,2 1 < $, where fi and (3 appear in (lty) . 

1 1 n 

(18) ^w\ losZ1= ^ E ^T^n^vw.---^). 

x 1 1 n>0 B lt - ,B n GC 1 — - 1 i=l 

oeuf =1 Bi 

which is analytic in a for \a\ < fi. 

Proof of Proposition^ We first show that for A C Z d , \a\ < /} and | | < 

(19) J2 \»l(B)\e m < 

sec(A) 

uniformly in A. Let S(B) be the image of B £ C(A), S'(.B) = {Ai, • • • , A;} c , where 
/ < Given Ai, • • • ,A„cA, such that U" =1 Aj is connected, we have 

1 7! n 

E ^ E« E ^ 1 7^n^ , (^)ii + ^i $1 ( A *) +$a ( A * 

seC(A) Z>1 ' fc x , •••,fe n n ' i=l 

S(B)={Ai,-,A n } c E?=ifci=« 



n oo 

n e tj^ii^)!! + ^ii $i ( a o + $2 (^ 
i=i fe=i 

n 

J~J( e Al|A'(A l )ll+^ll'f 1 (A l )+* 2 (A I )|| _ x y 



1=1 
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Therefore, 

n 

e a|a /4(5) = Y E e° E ^ l|Al1 ^(e^'MII^I* 1 ^* 3 ^)!! - 1), 

seC(A) n>l {A 1 ,.--,A„} C i=l 

BBx xgU" =1 AiCA 

uniformly in A C Z d . Together with Lemma [2] and translational invariance, this 
implies (1191) . Furthermore, it follows from ( fl5l) . ( |T9l) and translational invariance 
that 

(20) £ \^(B)\\x(B,B')\e & \ B \= ^ K(5)|e 5 ' B ' < |#|a, 

BeC(A) seC(A) 

uniformly in A G Z 6 *. 

Now, using (fl9l) and (|20|) . we inductively show that 

A n 

(21) i+$> e n^A(^)ii^---'^)i^ eai -'' 

n=2 B 2 ,- ,B„<=C{A) i=l 

for all iV > 2, uniformly in A C Z d . The case iV = 2 is trivially satisfied. Suppose 
that 

N-l n 

(22) i+j> £ ni^)n^---'^)i^ efii -' 

n=2 B 2 ,-,B„GC(A) i=l 

uniformly in A. It follows from (j!7p that 

N n 

Y, n e ni^(^)n^i,---,5n)i 

n=2 _B 2 ,-" ,B„ i=l 
N n 

=e e nK(^)i7^wiE n 

n=2B 2 ,-A«=l V ; ' GeC n (i,j)eE(G) 

Let G" be the graph one obtains after removing all the edges of G with endpoints 
at 1, and let (G\, ■ ■ ■ ,G n ) be a sequence of connected graphs of G' such that 
U k i=1 V(Gi) = {2, • ■ ■ , n} and n V(G f i ,-) = if i ^ j. Note that to each graph 

G' there corresponds fc! such sequences. Now, 

k 

i e n x(^-)i<e^i e nt n ^,^e n ^ 

GeG n (i,j)eE(G) k>l ' Gi,-,G k 1=1 (i,j)eE(Gi) G\ (i,j)e£(G|) 

where G\ runs over nonempty sets of edges with one endpoint at 1 and one in 
V(Gi). Furthermore, 

IE II X(B. 1 ,B,)\ = \ n (l + X (B u B t ))-l\< \x(BuB t )l 

G\ {i,3)eE{G\) ieV(Gi) ieV(G t ) 
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and hence 



N n 

1+ E n E \[\lA{Bi)MB x ,--- ,B n )\ 

n=2 B 2 ,-Ai=l 

k 

e nt n x(^) e ixc^x,^ 

fc>l ' Gi,-,G h 1=1 (i,j)eE(G t ) ieV(G0 



We now perform the sum over the graphs G\, • • • , by summing over partitions 
of {2, ■ ■ • , n} in sets V(Gi), % — 1, • • • , k, such that Yli=i — n — 1, and 



choosing connected graphs for each set of vertices. The number of such partitions 

{n-l)\ 

nil m, 



is J% . . It follows that 



N 



1+ E™ E HMWMBl'—iBn)] 

n=2 B 2 ,-,B n i=l 

TV 1 fc m; m; 

<1+Eu E II E \[\iA{B' i )\W{B' l ,-.^B' m )\Y,\ X {B l ,B^ 



k\ 

miH \-m k - 

mi mi 



n=2 m-,,- ,m k >l i = l B( ,-■ 4=1 1=1 



^ 1 + E^[ E n E liK(^ , )ii^i,--- ) K i )iEw^^) 

fe>l ■ l<mi,-,mfe<JV 1=1 B[,-,B' »=1 i=l 
1 AT m m 

= 1 + E^E E ni^)ii^,---,^)i^ix(5i,50i] fe 

fe>l ' m =l B[,-,B' m i=l i=l 



k\ 

k>l 



where we have used (120]) and (1221) in the last step. It follows by induction that 
(pjj) holds for all iV > 2. 

We now rewrite as an exponential. Consider sequences of connected graphs 
(G\, ■ ■ ■ , Gk) whose set of vertices form a partition of {!,••• , n}. Summing over 
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the number of vertices of each partition, and then over the partition, we get 

^X^ + EEf E rf^ x 

n>l k>l ' mx.-.mj. lli=i m Z- 
TiiH Ym^=n 

k mi 

x in e n^(A) e n xcs*,^)} 

i=l B lt - ,B mi i=l G£C mi (i,j)&E(G) 

1 fe mi 

= 1 + EE^ E II E n^w^-» B -«) 

n>l fc>l m 1 ,---,m fc >l / = 1 B\,--- ,B mi 8=1 

raiH \-nri} i =n 



1 + E^E E n^M^-.-^^r 

fc>l ' n>l B 1: ---,B n i=l 



Note that together with (JH]), it follows that |ZJ| < e |A| °. We now take the 
thermodynamic limit of well-defined quantities. Due to translational invariance, 

1 1 n 

in>g^-E E m^tII^m^---'^)! 

' ' n>0 b 1: - ,B„ec(A) ' 1—1 — i=l 

oeu? =1 Bi 

1 n 

' ' x£A n>0 B x ,— ,B„ i=l 

We now argue that the second term tends to zero as A / Z d in the van Hove 
sense; see for example [T7], Chapter 6. It follows from (fill and (12T1) that 

n 

E E ni^wii^-'-^j^a, 

n>0 B lr --,S n 8=1 

oeu™ =1 Bi 

uniformly in A 6 Z rf . For 77 > 0, 3 A' finite subset of 7L d such that 

n 

e e ni^)n^---'^)i<^ 

n>0 B 1 ,--- .B n i=l 

We let 

A(, := {x E A : A' + x C A}. 
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Then 

n 

TajEE E \[\^{B i )MB u ...,B n )\ 

' ' x£A n>0 B lt — ,B n i=l 

1 n 

= m(E+ E)E E ni^wii^'-'^)! 

' ' riGAf, zeAXAf, n>0 Bi,-,b„ i=l 

zeu^Bj^A 

|An| 1 

~ V \M |Aol 

Taking the limit r] — > and then A /* Z d , the RHS of the last inequality tends 
to zero, and it follows that 

1 1 - 

x 1 1 n>0B K ..,B„eC *- 1 — i=l 

which is analytic in a for \a\ < /L □ 

4.2. Proof of Theorem QJ 

Proof of (i). For A 6 Z d , we define 

1 TV Wa (e a| a Pa c -/3i^Hi (A)-/3 2 H fl2 (A) ) 

^A(a) := 7^J l0 S T -ftif« 1 (A)-/3 2 H fl2 (A) ' 

where v4 A appears in assumption (A5). For 5 = (A l5 • • • ,K\ B \) C e C, we define 



B 



a ^'(Aj), &(B) : = 
i=i j=i 



B 



= 1 = 2, 



and = $'(0) = 1. We have 



fc>0 ' XcA l>0 ' XcA 



-ciE e ^n*) n d+x(«.B j ))}x 

fc>0 n=l B lr -- ,s n eC(A) i=l l<*<i<« 

E?=i 1-8*1=* 

x £^e e (E H'^ l)! ri(- $1 ^)-^^)) n (i+x(5 m ,5n))}. 

£>0 ' j=l Bl- .Bj-eCfA) i=l l<m<n<j 
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Grouping terms with common support, and using (j!5p . one can rewrite the above 
expression as 

k \B„ 



e a|A|A Ae ~/3i 



n=0 ' Bi,— ,B n eC(A) m=l h>0 r=l s=l+k 



/3 2 $ 2 (A S ))] J] {l + x {B u B 3 )). 

l<i<j <a 



Therefore, 



F A (a) = i-{log^-logZr°}, 



where Z\ has been defined in ( 1161) . It follows from Proposition [T] that 

n n 

(23) limF A (a) = E £ -— — -(J]^)- J] ^(S,))^, • • • , SJ, 

x n>0 B 1 ,...,B„ — ~ i=1 i=1 

for |/?i )2 1 < /9 and \a\ < /i, such that lim A y Z d F A (a) is analytic in a for |a| < /U. 
We now define, for A C A', 

1 Tr H ( e a \^ e -^H Rl (K')-p 2 H R2 (K')\ 
Fa, A' (a) := r-rlog ' 



|A| Tr-u e ~ /3lHB i {A ' ) "' 32Hi? 2 (A ' ) 



For X, A C Z d , we introduce the function 

g(x,A)-- 



1, X C A 
0, otherwise 



and the measure on C(A') 
(24) 

1 |B| n |B| 

^' A,(jB) := <Wtt ^ ^nW^' A )) II (-A^(A S )-^2$ 2 (A S )), 

^ ^ n=0 r=l s=n+l 

for 5 e C(A'). We have 

1 " 
\F AA ,(a) - F A (a)\ < —(£ E II 1/^(^)11^1, ■ • • , ^01 + 



n>0 Bi,— ,B„eC(A') *=1 

u^^nA'VA^^ 



+ E E II^aTa^ii^i,---,^)!}. 



"■>0 B 1> — ,B„6C(A') i=l 
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Since A /* Z d in the van Hove sense, one can show using an argument similar to 
the one in the proof of Proposition [1] that, for \a\ < fl and < $, 

1 n 
1%^ iH E ^ II IfiU'WMBu ■■, B n )\ = 0, 

n>0 Sj,... ,s n eC(A') i=l 

and hence 
(25) 

1 n n 

^JT^'^E E H7r-Bj(II m)-II^)MBi,-, 5 n ) 

X X n>0 b 1 ,...,b„ 1 l ~ 1 — i=l i=l 

oeu? =1: B i 

which is analytic in a for \a\ < p,. 

To complete the proof of claim (i), we formally rewrite the moment generating 
function F using assumption (A5) as 

F(a)= lim lim ^-log^ (tt i (e a|ApA )) 

K/ldt^OO |A| 

= lim ^logw (e a|A| ^ +|A|i?A ) 
k/i d |A| 

1 Tr w ( e a Sx 6 A^'W+|A]flA e -Ef = ift^« i (A') 
= lim lim - — rlog( — = ) 



1 Tr w ,e a ^6A^W e -SLift^(A') 
= lim lim — — log( — — 5 ) + C, 

A/z<*AVz d |A| 6V Tr H ,e-S?=ift^(A') ; 

where the constant C appears in assumption (A5). Together with (125]) . this 
implies that, for \a\ < jl and < /3 



1 n n 

f h=e e i i , ? , R , i (ii ^) - n /a^m^, ■■■,B n )+ C, 

s 1j --.,s„ 



n>0 B 1 ,-,B n 1 — -' i=l i=l 



which is analytic in a. Claim (i) follows from analyticity of the moment generating 
function in a neighbourhood of the origin and Proposition 1 in [16j . 
Proof of (ii). 

The proof of claim (ii) is a straight forward extension of the result in [5] ; see also 
[TO] . We sketch the main steps of the proof. We rewrite the moment generating 
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function F as 



1 Tr H , e a E :c eA^+|A|«A e -Ei=iA^fl i (A') 
F(a) = lim lim — — ■ log( — —5 



where 



T •- 

A,A " Tr WA ,e a ^ e A^®i A , VA 



and 1a'\a is the identity on A'\A. We define the measure on C(A') 
for B G C(A'), and /Z AiA ,(</>) = 1. We have 

00 n 

(26) ^,a'=e^ e n^A'W n 

n=0 ' Bi,--,B Tl eC(A') i=l l<i<J<n 

Using the cluster expansion method, it is shown in [9] that, for | /3i ,2 1 < $ and 
a G ho s2 :={zeC: \$Sz\ < 



1 _ 1 ™ 
lim lim — logZ AA , = V V —Y\V a (B l )ip(B 1 , 

k/7L d tJ /TL d A 6 A ' A V ^ U?-1 11 
x /II n > Bl ,...,B n 1 »-l — 11 i=l 

0€Uf =1 Bi 



where ^ a (B) := liniA/'zdfimA' A*aa'- This implies the analyticity of F in a for 
a G 7 log 2 , which, together with the Ellis-Gartner theorem, implies claim (ii). 

□ 
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